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Abstract. Transformations achievable by linear optical components allow to generate the whole unitary group only 
when restricted to the one-photon subspace of a multimode Fock space. In this paper, we address the more general prob- 
lem of encoding quantum information by multiphoton states, and elaborating it via ancillary extensions, linear optical 
passive devices and photodetection. Our scheme stems in a natural way from the mathematical structures underlying 
the physics of linear optical passive devices. In particular, we analyze an economical procedure for mapping a fiducial 
2-photon 2-mode state into an arbitrary 2-photon 2-mode state using ancillary resources and linear optical passive Al- 
ports assisted by post-selection. We found that adding a single ancilla mode is enough to generate any desired target 
state. The effect of imperfect photodetection in post-selection is considered and a simple trade-off between success 
probability and fidelity is derived. 

PACS. 03. 67. -a Quantum information -03. 67. Lx Quantum computation - 42. 50.Dv Non classical states of the e.m. 
field, including entangled photon states; quantum state engineering and measurements 



9 L,: 
^— > ■ 

c ■ 
s : 



1 Introduction 

A quantum computer [ 1 1, although still a chimera as a concrete 
device, is already a venerable object for physicists, mathemati- 
cians and computer scientists, for the wide range of completely 
new perspectives that such a tool should offer for the develop- 
ment of science as well as for technological applications. 

Photon states are stable against decoherence, and are cur- 
rently produced and manipulated in modern laboratories. These 
features make the possibility of implementing quantum logic 
gates particularly attractive. One of the most promising archi- 
tectures for implementing a quantum computer by means of op- 
tical systems is based on a scheme proposed by Knill, Laflamme 
and Milburn (KLM) [2 1. In this scheme, information is encoded 
by (tensor products of) single-photon two-mode states of the 
quantized e.m. field; precisely, the qubit states are identified 
with a couple of single-photon states on two optical modes 
(dual rail logic) and multi-qubits are obtained by tensor prod- 
ucts. The basic ingredients for the elaboration of information 
in the KLM scheme are linear optical passive (LOP) compo- 
nents [3 1 — essentially, phase shifters and beam splitters — by 
which one is able to realize the single qubit gates; all other op- 
erations can be obtained in a non-deterministic way exploiting, 
in addition, ancillary optical modes and photodetection. One 
can show that, with the KLM scheme — hence, using only sin- 
gle photon sources, LOP devices and photodetectors — it is 
possible to simulate efficiently, i.e. by means of a polynomial 
amount of resources, an ideal quantum computer 1415 1 . 

It is worth noting, however, that in the KLM scheme, states 
that are not in the dual rail logic (e.g. the state |2000)) may be 



produced during the computation process, even if at the output 
they recombine to get back to the dual rail encoding. As it will 
be shown in the following, this is a consequence of the fact 
that the linear space spanned by all n-photon states (on a given 
number N of optical modes) is the carrier Hilbert space of an 
irreducible unitary representation of U(A^) which is associated 
in a natural way with the action of LOP devices. It seems then 
quite natural to investigate, in addition to the KLM dual rail 
logic, also the possibility of encoding information by means of 
rt-photon A^-mode states, with n > 1 and N > 2. The case 
where n = 1 and N > 2, with gates implemented only by 
LOP components, has been considered by Cerf et al. (6). This 
scheme is easily seen to be not scalable. 

In this paper, we will consider the case where information 
is encoded by n-photon states, with n > 1, on N > 2 modes 
and logic gates are obtained by LOP components and photode- 
tectors. As anticipated, this scheme stems in a natural way from 
the mathematical structures underlying the physics of LOP de- 
vices, structure that has been investigated in two previous pa- 
pers [7 8|. We will now address, as a first step, the following 
problem: to engineer any desired state — which may be re- 
garded as the 'input state' of a quantum computation process 
— in the chosen encoding space, starting with a fixed 'fiducial 
state', namely, a photon state that can be easily produced by 
single-photon sources. For the sake of definiteness, we will fo- 
cus on the case where n = N = 2. This is the simplest case that 
is not contemplated in the Cerf et al. and in the KLM schemes. 
Notice that in our case the building blocks of quantum infor- 
mation are qutrits instead of qubits. 
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The paper is structured as follows. In Section |3 the basic 
mathematical ingredients for a natural and systematic descrip- 
tion of LOP transformations are recalled. Next, in Sections |3] 
and our encoding and elaboration scheme is presented. The 
effect of realistic imperfect photodetection is then discussed in 
Section[5] Finally, in Section|6] we end up with some conclud- 
ing remarks. 



2 A group theoretical approach to LOP 
components 

A generic LOP transformation can be described as a 2 N— port, 
namely a black box with N input modes and N output modes. 
A pictorial representation is given in figure [J In two previous 
papers 17181 . it has been shown that the natural mathematical 
description of the action of LOP transformations is based on 
the theory of representations of semi-simple Lie groups and al- 
gebras; in this framework a special role is played by the Jordan- 
Sch winger map |9 10 1. In this section, we recall the basic in- 
gredients of such a description. Let us consider a set of N op- 



It is easy to show that the only maps with properties J3l4b are 
of the form: 



a, 
J 



b\ = M* 3 a] 



(5) 



where My is a N x N unitary matrix (M G XJ(N)). It is also a 
simple calculation to verify that a map of the form (0 preserves 
the canonical commutation relations: 



\bi,b 



,b j } = [bl,b]} 



SijI, 



0. 



(6) 
(7) 



Thus one can consider the two realizations of the (N di- 
mensional) Heisenberg-Weyl algebra given by {dj, a\, 1} and 
{bi, b\,l} and notice that by virtue of the Stone-von Neumann 
theorem 1 11 12 1 they are unitarily equivalent, that is, it exists 
an unitary operator U acting in the N— modes Fock space such 
that 



b] = WalU 



(8) 




Fig. 1. a black-box picture of a 2N— port based on L.O.P. transforma- 
tion 



Notice that the operator U is defined only up to an arbitrary 
phase factor. Since by construction U commutes with the total 
photon number operator this phase factor is fixed by the action 
of U on the vacuum state: 



U\0) 



|0). 



(9) 



This ambiguity can be removed if one considers an explicit 
construction of the unitary operator U. This can be done by 
means of the Jordan-Schwinger map. 



tical modes with the associated field operators 



i = l,2...N, 



(1) 



where the index i may label both spatial or polarization modes 
of the field, with the canonical commutation relations 



[aj,aj] = SijI, 



[ai,aj\ 



[a\, a]] 



0, 



(2) 



where I is identity operator. It is well known that the set of 
operators {a,, at I}, endowed with the canonical commuta- 
tion relations (|5J, are the generators of a realization of the 
N— dimensional Heisenberg-Weyl algebra W(N) |7j. We in- 
dicate with ftW) the bosonic Fock space associated with the 
chosen set of N modes. 

We are interested in Linear Optical Passive (LOP) transfor- 
mations, i.e., maps that are linear in the field amplitudes 



(3) 



(where the sum over repeated indices is assumed) and preserve 
the total photon number operator 



E 

= 1,..JV 



b] h = 



E 

i=l,...N 



a]di. 



(4) 



2.1 The Jordan-Schwinger map 

The Jordan-Schwinger (JS) map [9 H5), in its general formu- 
lation, maps a Lie algebra into an algebra of operators defined 
on a bosonic Fock space, this map being an algebra homomor- 
phism. The JS map is defined as follows. Let us consider a 
to— dimensional Lie algebra realized as an algebra of N x N 
matrices with a given basis of generators 



Q 



» = \\Q\f\\ a -1,2,... to i,j = 1,2, 
and commutation relations 



[Q( a >,QW]=c; 



,N (10) 



(ID 



Let us also consider a N— mode bosonic Fock space with field 
operators a«, aj i = 1, 2, . . . N and the (normal ordered) oper- 
ators 

dij = a\dj. (12) 
The operators d!2t satisfy the following commutation relations: 

[dij,dhk] = dikS^j — dhjSik- (13) 
One can consider the following set of bosonic operators: 



(14) 
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The map defined on the basis dlOi 

Q (Q) — -> JS(Q (a) ) 



{a) d- 



(15) 



extended by linearity, defines the JS map. It is easy to show 
that by virtue of the commutation relations Jl 31 the JS map is 
indeed an algebra homomorphism, namely 



[JS(QW),JS(QW)] = d°fJS(Q™). 



(16) 



Let us now come back to the transformation l|8}. The N x N 
unitary matrix Al can be written in terms of the exponential 
map as M = exp(A), where A is an element of the Lie algebra 
of the N— dimensional unitary group (namely a N x N anti- 
hermitian matrix). It is found that the related unitary operator 
U can be written in the following way by exploiting the JS and 
the exponential map, namely 

U = exp(J5(A)). (17) 

In order to check this, consider that, for e << 1, we have: 1 



U^cikU = exp{-eJS(A))a k exp(eJS(A)) 
~ a k + e[a k ,JS(A)] 



(18) 



where 

[ak,JS(A)} = Aij {a k a\aj - a\a 3 a k ^ 

(ia\a k + Sih)aj ~ ajajakj = A k3 a r 



(19) 



The JS map allows to fix the arbitrary phase factor in (|9j. In 
fact, since the JS(A) is a normally ordered operator, we have: 



exp(J5(A))|0) = |0), 



(20) 



so that = 1. 

To summarize we have shown that LOP transformations on 
TV modes are described by means of the TV— dimensional uni- 
tary group acting on field operators as in (|5j- Such an action of 
the TV— dimensional unitary group induces a bosonic represen- 
tation of the group U(iV) acting on the N modes bosonic Fock 
space 

U = T {N \M) (21) 

that can be explicitly defined by means of the JS map. Indeed 
it is easy to check that 



rW(MiM 2 ) = T {n) {M 1 )T {n \M 2 ). 



(22) 



in correspondence with the decomposition 



(24) 



n— 0.1, ...oo 



where Tin is the subspace with n photons on N optical modes. 
A simple calculation shows that the following relation holds: 



dimftW = 



Hence, the subspace Tin can be seen as the space of a qu-dit 
with d = dirriH„ . The (sub) representation with n = is the 
trivial representation of V(N): 



(n + N-l)\ 
n\(N-l)\ 



(25) 



r o (Ar) (M)|0)H0}; 



(26) 



while a special role is played by the n — 1 (sub) representation 



since 



(0|o fc l7o||0) = (0\U(U*a k U)a}\0) = (0\b k al\0) 

= J2 M kh{0\a h aJ\0)=M k i, (27) 



where we have used the fact that UW = I and W\0) = 

U\0) = |0). The nature of the representations {ri JV ' ) }n=o,i,... 
has been studied in detail in 17181 . A remarkable result is that 
each (sub) representation Yn is a irreducible unitary repre- 
sentation (IUR) of the group U(JV). For N = 2 and n = 1, 

(2) 

Y± is the relevant (sub) representation for the implementation 
of single qubit gates in the framework of dual rail logic 112 17 1 . 



3 About two-mode multiphoton states 

For the sake of definiteness, in the following we will focus on 
the case where N ~ 2. This configuration is at the basis of the 
KLM scheme for a quantum computer |2), in which the qubit 

Hilbert space is identified with the space n\ ' of one photon 
on two modes. An important requirement for a well defined 
quantum computation is the ability to perform an arbitrary one 
qubit gate [ 13 1, that is, a generic unitary transformation in the 

qubit space Hi . It is a well known result that in KLM scheme 
every one qubit gates can be implemented with only two-modes 
LOP transformations. This follows directly from the fact that 

Tf is the fundamental representation of the group U(2) acting 

(2) (2) 

on the one photon subspace Ti\ (with dimTi\ = 2). This 
is no more true in those subspaces characterized by a larger 



Since, by construction, Y^ N \M) commutes with the total 
photon number operator, the unitary representation yW can number of photons. For n > 1, Tn' is a spin— ^ representation 
be written as the direct sum of unitary (sub) representations 
acting on the subspaces with fixed photon number 



(23) 



n— 0,1, ...oo 



,(2) 
n 

(2) j— ■ ~ (2) 

acting in the n photon subspace TL n |7| (with dim7i;, = 
n + 1). Thus, in the case where n > 2 it is no more possible 
to realize a generic unitary gate in the n— photon subspace and, 
in general, it could not exist a LOP transformation (associated 
with some unitary matrix M) such that 



1 For a rigorous proof one can use the well known formula 
e x Ye~ x — exp(adx)^, for linear operators X and Y. 



\<P)=YW(M)\iP ) 



(28) 
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for a generic couple of normalized state vectors \i/jq), \4>) € 

(2) 

Tin ; in other words, for n > 2 not all the normalized vectors 
belong to the same U(2)— orbit. 

Let us recall that, given a representation T of a group G in 
a Hilbert space H, the orbit O^ of the group passing through 
a given vector \ipo) E H is defined as the set of all vectors 
IV*) € T~i sucn that = Y(g)\ipo), for some g 6 G. In the 
case where n = 1, the orbit of the group U(2) in passing 
through a vector of unit norm fulfills the whole unit sphere in 

(2) 

H\ .In the multiphoton case, the orbit of the group U(2), act- 

(2) (2) 

ing in Tin via the representation T„ , passing through a nor- 
malized state vector \ipo) E Tin , is only a proper sub-manifold 
of the unit sphere. 

In order to illustrate these arguments explicitly, let us con- 
sider a generic SU(2) matrix: 



M = 



a 
-(3* a* 



(29) 



with a = 



cos ( 



(3 = e 1 ^ sin 9. The two-photon sub- 



space T-6p has dimension d = 3, hence it can be seen as a 
qutrit space. In the number basis {|20), |11), |02) } the operator 

(2) 

T 2 (M) has a matrix representation 



r 2 (2) (A/) 



a 



V2a*(3 



-V2a*f3* a* 



(30) 



It explicitly shows that it is not possible to realize every (qutrit) 
unitary transformations. Also notice that 

~(2) 



r^(M)\ll) = V2a/3|20) + (H 2 - 



) |ll)-V2a*/r|02>, 
(31) 

from which it is apparent that the state vectors 1 1 1) and 1 20) (or 
1 02)) do not belong to the same orbit. 

To summarize, in the multiphoton case two categories of 
problems arise that are not present in the single photon case: 1) 
given a state vector \tpo), there is in general no LOP transforma- 
tion that allows to map \ipo) into an arbitrary target state vector 
\ip); 2) it is not possible to perform every qutrit unitary gate 
only with two-mode LOP transformations. The latter problem 
was investigated from different points of view in 1 1 41 1 5 1 1 61 . 
These problems are related to the DiVincenzo's criteria 1131 
for a well defined quantum computation, namely the point 1) is 
related to the ability to initialize the state of the qutrit to a sim- 
ple fiducial state; the point 2) is related to the ability to perform 
a universal set of quantum gates. In the following sections, we 
consider the first problem and suggest a solution based on pho- 
todetection on ancillary modes and conditional post-selection. 



4 Projection via a post-selection protocol 

A remarkable property of IURs is that every orbit is total 1171 . 
This means that given a normalized target state vector \<j>) £ 

Tin and an orbit O^ of the IUR T„ , it is always possible 
to find a e Otf, with a non vanishing projection along \<f>): 

W#> * 0. (32) 



While the existence of a non-vanishing projection follows from 
the properties of the irreducible representations T„ \ how to 
realize physically (at least in principle) such a projection is a 
matter of a different nature. In the following we discuss with 
some examples a procedure based on photodetection on ancil- 
lary optical modes and conditional post-selection. The result is 
a non-deterministic protocol that allows to map a fixed input 
state into a desired target state with a certain probability. In 
order to illustrate the idea, let us consider the case of a qutrit 

(2) 

encoded in the subspace with two photons on two modes ri 2 . 
The proposed procedure consists in four steps. The first step is 
to initialize the qutrit system in a fixed input state |^> )- The 
second step is to add one extra optical mode that plays the role 
of an ancilla: the state of the ancillary mode is initialized in 
a number state with m photons and we consider the extended 
state 

|V>o) — > W)\m). (33) 

(3) 

Hence, the relevant space for the system+ ancilla is n^+m - The 
third step is to perform a three-mode LOP transformation that 



acts on H^K 



via the IUR T. 



(3) 



|^o)|m)— ftr|^ o )|m) = |^ o )|O) + |0 1 )|l) 



• \m). 
(34) 
the 



The final step is a post-selection on the ancillary mode: 
target state is obtained, with a certain probability V m , in corre- 
spondence of the detection of m photons on the ancillary mode. 

Overall the transformation of the initial state is described 
by a completely positive map £ ( m > which depends on the initial 
preparation of the ancilla mode. The Kraus-Sudarshan form of 
— > £( m ) is of course given by 

IVoXVol — > e^\\^){M) = E ^VoWol^ - 

m' 

(35) 

where — (m'\U\m). The post-selection conditioned on 

the photodetection of m' photons on the ancillary mode corre- 
sponds to a single branch of the map i.e. to the transformation 



hM — A% ) \i> ) = \ct> m ,}. 



(36) 



In the following two examples are presented with m = 0, 1. In 
the appendix|X|it was shown that adding one ancillary mode is 
indeed sufficient in order to obtain the optimal working point. 



4.1 On the ability to initialize the state of a qutrit to a 
simple fiducial state 

Let us consider two computational modes with one extra ancil- 
lary mode and a three-mode LOP transformation: 



Mt. 



,j = 1,2,3. 



(37) 



Let us also take the third (ancillary) mode initialized in the 
vacuum state. Following the procedure outlined above, here 
we answer the question of whether is possible to find a three 
modes LOP transformation such that, after a photodetection on 
the third (ancillary) mode, a generic qutrit state 



A|20) +B\11) +C\02) 



(38) 
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Fig. 2. Photodetection assisted three-mode L.O.P. transformation 

is obtained, with a certain probability, on the first and sec- 
ond (computational) modes. As input state we select the state 
|-0o) 12 = | Il)i2 that will be extended with one ancillary mode 
initialized in the vacuum state (to = 0) 



|ll)i2|0) a 



(39) 



The subscripts indicate the mode labels and will be omitted in 
what follows. 

The action of a LOP transformation acting on j39t yields 

to: 

\f) = M p iM 32 <44|0Q0) p,q = 1,2,3. (40) 

The global (three-mode) output state, obtained after the three- 
mode LOP transformation has the form: 



(41) 



|V>) = |0o>|O> + |0i>|l> + |02>|2>, 

where \<j> n ) are two-mode states. A post-selection conditioned 
to the vacuum on the third optical mode gives: 

|V> — > |0o)|O) =M p iM, 2 at a t|00)|0) p,q= 1,2 , (42) 
where 

|0o) = n/2MuMi 2 |20) 

+ (M n M 22 + M 2 i Afu) |H) 

+ V2M 2 iM 22 |02) (43) 

is the un-normalized two-mode output. The square modulus 
"Po = (0o 1 00 ) gives the probability of success of the vacuum 
measurement. From a mathematical point of view, the question 
is whether is possible to find, for every target state J38l >. an 
unitary matrix M such that the output state d43t and the target 
state J38l > are equal apart of a normalization (and phase) factor. 
The following propositions hold: 

Proposition 1 for any a, (3, 7, 5 with \a\ 2 + 
exists an unitary matrix 



< 1, there 



M 



a 7/ A e 3 
(3 S/k e 4 
ei e 2 /k e 5 



(44) 



for some e\ , . . . e§ and real k =/= 0. 

Proof: in order the matrix M to be unitary the following equa- 
tions have to be satisfied: 



M 2 - 




"leil 2 


= 1 


(45) 


l7| 2 


f l<5IN 


-hi 2 


= k 2 


(46) 


leal 2 4 


|e4|N 


-N 2 


= 1 


(47) 


* 

a 7 


f /?*<$- 


- e^e 2 


= 


(48) 


a e 3 4 




- ejes 


= 


(49) 


7 e 3 ^ 


- 5*e A - 


- e 2 e 5 


= 0. 


(50) 



Let us suppose that e\e 2 7^ 0. Equation (145 \ and J48l > yield to 

\a*-f + (3*6\ 



e 2 



that inserted in (I46> gives 

fc 2 = l7| 2 + H 2 



|a*7 + /3*^| 2 
l/3| : 



1 



(51) 



(52) 



Once k, e\ and e 2 are found, the remaining coefficients can be 
easy computed by an orthonormalization algorithm. Otherwise, 
if e\e 2 = 0, there is always the trivial solution e\ = e 2 = e 3 = 
e4 — and e5 = 1. □ Notice that one can always choose 
a, j3, 7 and 5 such that the following normalization condition 
holds: 



2|o7| : 



*<5 + /3 7 | 2 + 2|/3^ 



1. 



(53) 



The previous proposition implies that starting from the state 
1 11) 1 0) for any normalized target state (1381 there exists a LOP 
three-mode transformation such that, after a post-selection mea- 
surement corresponding to the vacuum on the ancillary mode, 
the following transformation is obtained 



111) 



(54) 



Proposition 2 Vq — k 2 is the probability of success of the 
post-selection measurement 

Proof: within the scheme of figure[2]the output state is 



|0 O ) = - (\72a 7 |20) + (aS + (3j) 1 11) + V2(3S\02) 



(55) 



With the normalization condition d53l we obtain that Vq = 
(0o 1 0o ) = fc~ 2 is the probability of success. □ Given a nor- 
malized state vector in the form (1381 the optimal gate corre- 
sponds to the maximum of To (or the minimum of k 2 ) with 
constraints: 



V2q!7 = A 
aS + (3j = B 
V2/3S =C 



(56) 



4.1.1 Examples 



Let us suppose that we want to reach the state 1 20) starting from 
|11): 

|11) — |20). (57) 

In the following we are going to describe in which way the 
transformation d57i can be obtained with the maximum prob- 
ability. Notice that the state 1 20) is obtained from d38l > taking 
B = C = 0, thus we are now looking for U(3) matrices of the 
form 



M 



a 7/fe e 3 
e 4 
ei e 2 /k e 5 



(58) 
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The normalization condition (I53> implies that 2 1 0:7! 2 
d52t reads as follows 



+ 



I - a 



1, and 



(59) 



Hence, the maximum of probability is P mnl = 1/2 (that corre- 
sponds to the minimum of k 2 ) and it is reached for \a\ 2 = 1/2. 
Notice that this is the maximal probability allowed in the given 
set up 1 18 1. The corresponding unitary matrix can be chosen as 
follows: 



M 



.1 L 

V2 V2 

1 

_J L 



(60) 



This is not the only solution, with this choice the three-mode 
gate (I60> can be decomposed as product of two-mode gates in 
the following way 



M 



10 
1 
1 



1 L n 

-1 L n 

1 



(61) 



The circuital implementation is schematically represented in 
figure [3] and consists of a symmetric 50% beam splitter on the 
first and secondmode (9 — 7r/4, = 7r/2), followed by a swap 
operation between the second and third mode. 



8=7t/4 
*=Jt/2 



b, 



Fig. 3. Circuital implementation scheme of the three-mode L.O.P. 
transformation 1601 with post-selection procedure 



As an other example we are going to describe a post— se- 
lection assisted LOP transformation with one ancillary mode 
which is initialized with one photon. The computational space 



is H.^' with an ancillary space TOf 1 , hence the global space is 

(3) 

Ti. 3 . With a procedure analogous to that presented above, it is 
easy to shown that the same circuit of equation ( 1601 (and figure 
13 allows to perform the transformation 

|20>|1> — 111)11), (62) 
with an optimal probability of 1 /2. 

5 Effects of imperfect photodetection 

In the previous sections we have made the assumption that 
all the components are ideal: in this section we discuss the 
presence of real photodetectors. Light is revealed by exploit- 
ing its interaction with atoms/molecules or electrons in a solid: 
each photon ionizes a single atom or promotes an electron to 



(i) 



a conduction band, and the resulting charge is then amplified 
to produce a measurable pulse. In practice, however, available 
photodetectors are not ideally counting all photons, and their 
performances are limited by a non-unit quantum efficiency 77, 
namely only a fraction £ of the incoming photons lead to an 
electric signal, and ultimately to a count. For intense beam of 
light the resulting current is anyway proportional to the incom- 
ing photon flux and thus we have a linear detector. On the 
other hand, detectors operating at very low intensities resort 
to avalanche process in order to transform a single ionization 
event into a recordable pulse. This implies that one cannot dis- 
criminate between a single photon or many photons as the out- 
comes from such detectors are either a click, corresponding to 
any number of photons, or nothing which means that no pho- 
tons have been revealed. These Geiger-like detectors are of- 
ten referred to as on/off detectors. For unit quantum efficiency, 
the action of an on/off detector is described by the two-value 
POVM {n Q = |0) (0|, 77i = I - 7T }, which represents a par- 
tition of the Hilbert space of the signal. In the realistic case, 
when an incoming photon is not detected with unit probability, 
the POVM is given by fEfl 



n (v) = J2(i-v) k \k){k\, 



fe=0 



tfifa) = i-n ( v ), 



(63) 



with 7/ denoting quantum efficiency. As a consequence the con- 
ditional state, occurring when the event "no click" is registered 
is no longer the pure state given in Eq. ( I43> . The conditional 
state is now given by the mixed state 

go = — Tr 3 [I^X^l I <8> I <8> JT (»?)] 

-MJ 



-Ed 

Po ^ 

' k=0 



(64) 



where \ip) is given in Eq. ( 1401 . \<j>k) are the unnormalized states 
corresponding to an ideal (unit quantum efficiency, perfect dis- 
crimination) photodetection of k photons and Po is the global 
probability of the "no click" event, i.e 



Pn 



2 

£ 

k=0 



(65) 



The (unnormalized) conditional state |</) ) is given in Eq. d43t 
whereas \<f>k), k = 1, 2 are given by 

\(fn) = (M n M 32 + MaiMia) |10) 

+ (M31M22 + M21M32) |01), (66) 
|0 2 ) = V2Af 3 iM 32 |00) . (67) 

Realistic photodetection thus degrades the quality of the prepa- 
ration. In order to asses the whole procedure we use fidelity to 
the target state i.e. 



F 



1 



(0o I ) 
1 



(0o|0o)Po 



(0o I go I <Po) 

x;i<0oi0*)i a (i -»?)*. 



(68) 
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Since the conditional states \cj>k) 
obtain 

(00 l^o ! 



are mutually orthogonal we 



F 



Vo 
Pa' 



(69) 



Therefore there is a simple trade-off between the probability of 
success and the quality of the preparation, which can be used 
to suitably adapt the procedure to the desired task. 

In the case of postselection corresponding to a click of the 
photodetector the roles of 77o and 27i in (I63> are inverted. A 
click on the ancillary mode corresponds to the preparation of 
the computational modes in the mixed states 



where 
and 



A k ( V ) = l-(l-r,) k 



k ■ 



(70) 
(71) 
(72) 



k=l 



The corresponding fidelity to the target state is 



F = 



(0x10!) (4>MPi 

which simplifies to 

Ai(r])(ct>i\<t>i 



Y^MvWMl 2 , (73) 



F = 



77 A ' 



(74) 



Hence, also in this second example a simple trade off between 
probability of success and fidelity of real processes is obtained. 

In general, the probability of success and fidelity are inde- 
pendent quantities in the sense that the maximization of the 
success probability does not imply the fidelity optimization. 
For example, the optical circuit in figure [5] corresponds to the 
maximal probability of success for both the transformations 
|11)|0) -> |20)|0) and |20)|1) -> |11)|1) with an optimal fi- 
delity for the former and a non-optimal fidelity for the latter. 



6 Conclusive remarks 

In this paper we have addressed the problem of whether in addi- 
tion to the KLM dual-rail quantum computation one can con- 
sider a more general n-photon iV-mode encoding scheme; in 
other words, whether there is room for quantum information 
processing based on multiphoton encoding of qudits. In partic- 
ular, we investigated the problem of the system initialization 
in Hilbert spaces that are carrier spaces of irreducible unitary 
representations of unitary groups, representations which are as- 
sociated in a natural way with LOP transformations. Focusing 
on the case of the 2-photon 2-mode encoding, we found that 
LOP devices assisted by post-selection measurements allow to 
engineer any desired state in the encoding space starting from a 
suitable fiducial state; moreover, we have shown that the use of 
a single ancilla mode is enough to ensure the maximum prob- 
ability of success. The effects of imperfect photodetection in 



post-selection have been considered and a simple trade-off be- 
tween success probability and fidelity has been derived. 

Of course the lack of further generality and detail in our 
present investigation is something to be remedied in the future. 
However, we think that it would unrealistic and may be futile, 
at this preliminary stage, to try to solve in its full generality 
the problem of simulating an ideal quantum computer within 
the encoding scheme that we have proposed here. Our main 
purpose is to suggest that a deeper understanding of the mathe- 
matical structures underlying LOP devices could be a powerful 
tool for the further development of optical quantum computa- 
tion. 



Acknowledgments 

We wish to thank Prof. G. Marmo of the University of Napoli 
'Federico IF for his invaluable scientific and human support. 
This work has been supported by MIUR through the project 
PRIN-2005024254-002. 



A One ancilla mode is enough 

In the body of the paper we analyzed in some details the prepa- 
ration scheme based on a single ancillary mode. In this section 
we show that adding a single ancilla is enough in the sense that 
with multiple ancillary modes no improvements of the proba- 
bility of success can be reached. We consider the case in which 
the initial input state is the two photon state |11) and discuss 
the m ancillary modes generalization of the propositions[Oand 
[2] The matrix (I44l > has the following generalized expression in 
the case of m ancillary modes: 



M 



j/k 
6/k 



e 3 
e 4 



ei J e 2 I /kE 5 



(75) 



where ei are m-component complex vectors and E§ is a m x m 
matrix. Equations ( 145 \ J46t and J48t become: 



a*a + [3* [3+ |ei | 2 = 1 
7* 7 + r«5+|e 2 | 2 =k 2 
a*~/ + (3*5 + (ei,e 2 > = 0. 

Taking (ei, e 2 ) 7^ Owe obtain 

l"*7 + /3*<5| 2 



P = | 7 | 2 + H 2 + 



where 



ei,e 2 



cos0|(l - |a| 



|ei||e 2 | cost 



(76) 
(77) 
(78) 



(79) 



(80) 



From i79i it follows that the maximum probability is reached 
at I cos 9\ — 1 and correspond to the value in i5H . Otherwise, 
in the case (ei,e 2 ) = 0, there is always the trivial solution 
ei = e 2 = e3 = e 4 = and E$ = I. 
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